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Studying the problem of wave propagation in media with resistive boundaries can be made 
by searching for "resonance modes" or free oscillations regimes. In the present article, a 
simple case is investigated, which allows one to enlighten the respective interest of different, 
classical methods, some of them being rather delicate. This case is the ID propagation 
in a homogeneous medium having two purely resistive terminations, the calculation of the 
Green function being done without any approximation using three methods. The first one is 
the straightforward use of the closed-form solution in the frequency domain and the residue 
calculus. Then the method of separation of variables (space and time) leads to a solution 
depending on the initial conditions. The question of the orthogonality and completeness of 
the complex- valued resonance modes is investigated, leading to the expression of a particular 
scalar product. The last method is the expansion in biorthogonal modes in the frequency 
domain, the modes having eigenfrequencies depending on the frequency. Results of the three 
methods generalize or/and correct some results already existing in the literature, and exhibit 
the particular difficulty of the treatment of the constant mode. 

PACS numbers: 43.20Ks, 43.40Cw, 02.30Tb, 02.30Jr 



I. INTRODUCTION 

Studying the problem of wave propagation in me- 
dia with resistive boundaries can be made by search- 
ing for "resonance modes" (see Ref.l), or free oscilla- 
tions regimes. These modes can be non-orthogonal for 
the ordinary scalar product, entailing some difficulties 
depending on the mathematical treatment, made either 
in the time or frequency domain. Two classical methods 
exist for such a problem, and can be used either for a 
scalar, second order differential equation, or for a system 
of two equations of the first order. They have been es- 
pecially used for the problem of a ID medium with one 
resistive boundary, the other boundary condition being 
of Dirichlet type: 

i) in the time domain, the use of time and space vari- 
able as separate variables leads directly to the basis of 
modes, but they are non-orthogonal for the most com- 
mon product, and difficulties occur when searching for 
the coefficients depending for instance on initial condi- 
tions. Nevertheless, for a particular case, Oliveto and 
Santini, 2 and Guyader 3 have solved the problem, and 
Rideau, 4 using a system of equations of the first order, 
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found a scalar product making the modes orthogonal (see 
also Refs.5-7), and gave the proof of completeness. 

ii) in the frequency domain, the equations to be solved 
are ordinary differential equations with boundary condi- 
tions depending on frequency, but the use of orthogonal 
decomposition is possible. This leads to eigenmodes and 
eigenfrequencies depending on frequency. It is the case for 
the classical theory of room acoustics (see e.g. Morse and 
Ingard 8 ), using biorthogonality. To return to time do- 
main in order to deduce the resonance modes is a rather 
delicate task, especially because of the calculation of the 
derivation of eigenfrequencies with respect to frequency. 
Biorthogonality has been used also for duct modes (see 
e.g. Ref.9) Another approach has been recently used by 
Trautmann and Rabenstein, 10,11 using a system of first 
order equations (these authors treat the case of two re- 
sistive boundary conditions). 

The present article is devoted to the study of the sim- 
ple ID case, when the two boundaries are resistive. One 
goal is to exhibit how the different methods articulate. 
We start by using the fact that a straightforward solution 
exists for the wave equation with source, by applying the 
residue calculus to the closed-form of the Fourier domain 
solution: as discussed by Levine, 12 this closed-form solu- 
tion, avoiding the sum of a series, is "relatively poorly, 
if not entirely, unknown to the general acoustics com- 
munity" . All calculations can be carried out analytically 
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without any approximation, exhibiting the properties of 
the different methods (however many previous papers re- 
strict their content to small impedance, or admittance, 
at one extremity, using perturbation methods). The case 
under study corresponds to one-dimensional propagation 
in a homogeneous medium bounded by two other semi- 
infinite media with different characteristic impedances, 
dissipation being therefore due to radiation at infinity, 
ft is especially interesting because of its physical signif- 
icance (it is probably the simplest radiation problem), 
and also because it realizes one of the possible transi- 
tions between Neumann and Dirichlet boundary condi- 
tions. Notice that in the context of optics and quantum 
mechanics, the problem has been studied including the 
outside media by Leung et al, 13,14 the resonance modes 
being called quasinormal modes. 

In section II, the equations to be solved are stated, 
with some possible physical interpretations. As a first 
step, the classical, closed-form solution of the Green 
function in the frequency domain is established (section 
III), with its inverse Fourier Transform, corresponding to 
the successive reflections (section IV). The second step 
is the residue calculus in order to determine the reso- 
nance modes (section V, the basic result being given by 
Eqs. (26)). Then results are compared to those of the 
two aforementioned methods, i.e.: i) the method of sepa- 
ration of variables (section VI), which gives the result for 
given initial conditions (the corresponding results being 
Eqs. (37), (48) and (49)); in this section, the question 
of orthogonality and completeness of the modes is inves- 
tigated, ii) the method of eigenmodes in the frequency 
domain (section VII). For the two methods, both sec- 
ond order scalar equation and first order system of two 
equations are used successively, with emphasis on the ex- 
istence of a constant mode. 

II. STATEMENT OF THE PROBLEM, PHYSICAL 
INTERPRETATION 

The Green function g(x,t | x 0} t ) for the wave equa- 
tion is solution of the following equation: 



dt t ]g(x,t) = -5(x-x Q )8(t-to) (1) 



where x and xo are the spatial coordinates of the re- 
ceiver and source, respectively (or vice- versa), t and to 
the times of observation and excitation, respectively, c 
the speed of sound. 5(x) is the Dirac function. 

For sake of simplicity, xo and to are considered to be 
fixed. Moreover in the whole paper, the choice of to = 
is made. For negative t, the function is zero, as well 
as its first derivative. The Green function satisfies the 
following boundary conditions: 

c(d x g(x, t) = d t g(x, t) at x = 0, (2) 
cQd x g(x, t) = -d t g(x, t) at x = I (3) 

where £ = ZJ pc, p is the density of the fluid, and 
Z the impedance at x = 0, which is assumed to be a 
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FIG. 1. ID medium bounded with two other media. 



real quantity, independent of the frequency. Similarly, 
Q = Zi/pc , where Zi is the impedance at x = I (I being 
positive) . 

An obvious physical interpretation for quantities ( and 
Q is the following: consider for x < and x > £ (see 
figure 1) two media with characteristic impedances p_c_ 
and p+c+, respectively. If the media are non dissipative, 
impedances are real, and can be larger or smaller than 
the impedance of the bounded medium, pc. Moreover, 
they are positive, because they correspond to waves out- 
going from the bounded medium. Therefore this is the 
problem of planar pressure waves in a stratified medium, 
the direction of propagation being normal to the inter- 
faces. A generalization to more complex stratified media 
would be possible, at least numerically. In this problem, 
the Green function corresponds to the acoustic pressure: 
of course, it has not the dimension of a pressure, but 
the solution for a "concrete" problem with source can 
be easily solved, as explained in standard textbooks, and 
discussed in a recent paper by Levine. 12 

Other problems correspond to the previous equations: 

i) in an approximate way, ignoring higher order duct 
modes, the problem of planar waves in a rigid walled duct 
terminating in two semi-infinite ducts with different cross 
sections areas, the quantities ( and Q being the ratios of 
the areas. The approximation is good at low frequencies. 

ii) the problem of a dissipative termination : the ter- 
minal impedances Z and Zi can be the impedances of 
dissipative media (at low frequencies, a porous medium 
open to a large space can be an approximation of a pure 
resistance, due to viscous effects). 

In all the previous problems, the quantities £ and (e 
are real and positive, the terminations being passive. For 
active terminations, they can be negative. An example is 
the beginning of self-sustained oscillations in musical in- 
struments: a nonlinear excitator, like a reed for a clarinet, 
can be linearized as a pure resistance. When the main 
control parameter, i.e. the pressure in the mouth of the 
musician, increases, the resistance becomes negative, the 
static regime becomes unstable, and an oscillation starts 
as an increasing exponential (see e.g. Refs.f 5,16). 
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Obviously analogous problems for mechanical vibra- 
tions or other wave fields are numerous. 

III. CLOSED-FORM SOLUTION FOR THE FOURIER 
TRANSFORM 

The Fourier Transform (FT) of g(x, t) is denoted 
G(x, u>) (throughout the article, functions of time are 
written in small characters, and their FT are written in 
capital characters). It is equal to: 

/+oo 
g(x,t)e- luJt dt, where (4) 
-oo 

i r + °° 

9{x ' t)= 2^J G{x,uj)e wt duj. (5) 
The FT of Eqs. (1) to (3) arc found to be: 

{d 2 xx +w 2 /c 2 )G(x,cj) = -6(x-x ), (6) 
cC,d x G(x,uS) = iujG(x,uj) at x — 0; (7) 
cQd x G(x, u>) = —iuG(x, u>) at x — £. (8) 

While terminal impedances are independent of fre- 
quency, boundary conditions are frequency dependent. 
Nevertheless a classical, closed-form, solution is already 
known, which has been especially used in Ref.17. If 
x ^ x solutions of Eq. (6) can be written as: 

G(x, to) = A~ cosh [iujx/c + rj\ ifx<x ; 
G(x, u) = A + cosh [iuo{£ — x)/c + rje] if x > x n . 

For the boundary conditions, the following definitions 
are used: 



( = coth?7 ; r = e 



- e -tv - 



coth?7£ ; ri = e 



(C-i)/(C + i) 

(0-l)/(0 + l), (9) 

where r and are the reflection coefficients. The quan- 
tity r] satisfies: 2rj = — (In \r\ + i arg(r)) [2tt]. Because r 
is real, we choose the following definition: 



V = Vr + ; ji = or 1. 



(10) 



Two cases exist: i) if \(\ > 1, r > 0, fi = 0; ii) if \(\ < 
1, r < 0, ii = 1. Similar remark and definition can be 
applied to boundary x = I: 

m = ntr + «mW 2 ; W = or 1. (11) 

The case £ = 1 (semi-infinite tube or medium) corre- 
sponds to r\ = oo: it is discussed in the next sections. 
Except for the last one, most of the following calculations 
are valid for all cases. At x = x , writing the continuity 
of the function and the jump of its first derivative, the 
following result is obtained: 

c cosh [rj + iujxo/c] cosh [% + iu>(£ — x)/c] 



G(x, 



iw 



smh(iLu£/c + r] + 



(12) 

if x > x and a similar result if x < x , by interchang- 
ing x and xq. 



IV. SOLUTION IN THE TIME DOMAIN 
(SUCCESSIVE REFLECTIONS) 

Eq. (12) can be transformed in the time domain, lead- 
ing to a solution corresponding to the successive reflec- 
tions of the Green function in infinite space at the two 
boundaries. It will be the reference solution for the check 
of the validity of the modal expansion. The sinh function 
of the denominator can be written as 



smh(iuj£/c + r] + t]e) = 



^ _ e -2)j-2ijf-2iw£/c 



and, if the modulus of the exponential at the denomi- 
nator is less than unity (this is discussed hereafter), as: 

sinh-^t^/c + f] + vt) = 2e- r >- r 't-™e/c 

[1 + F(uj) + F 2 (cu)+F 3 (uj) + ...]. (13) 

F(u) = exp(— 2?7 — 2rj£ — 2iuj£/c) is the function corre- 
sponding to a complete round trip of a wave in the tube, 
of duration 2£/c. Concerning the numerator of (12), it 
can be written: exp(+r] + r]i + iLLi£/c) G p (x,uj) c/4, where: 



G p (x,oj) = e -M*-*o)/c + re -M*+*o)/c 
+r£e~ iui * 2l ~ x ~ xa ^ c + rrie~ iui * 2l ~ x+xa ^ c . (14) 
Therefore the Green function is: 

G(x, lu) = ^-G p (x, u) [1 + + F 2 (lu) + ...]. (15) 

The factor G p (x,w)/iw corresponds to the four "pri- 
mary" waves arriving during the first cycle of dura- 
tion 2£/c, and this packet is simply reproduced at times 
2£/c, 4£/c, Q£/c, etc... (see for a detailed explanation 
e.g. Kergomard 15 ). The inverse FT of the function 
Gp(x,uj)/iu>, denoted h p (x,t), is obtained by taking into 
account the zero condition for negative times. The result 
is found to be, whatever the sign of (x — xq): 

h p (x,i) =H[t-\x-x \ /c] +rH[t- (x + x )/c] + 
rH [t-{2£-x- x )/c] + rr e H [t - (2£ - \x - x \)/({\16) 



where H(t) is the step function. Finally 

f (t) + (f * j 

f(t)=rr e 8(t-2£/c). (18) 



g(x, t) = -h p (x, t) * [S(t) + f{t) + (/ * f)(t) + ...] ; (17) 



Condition of validity of expansion (13) is \m\ < 1. 
We notice that if ( is real and positive, \r\ < 1, and 
similarly for Q. Therefore the condition is satisfied when 
the two boundaries are dissipative, or, more precisely, if 
the combination of the two reflections is dissipative. The 
case \rre\ > 1 will be discussed in section VC. Other 
comments can be made: 

- the article is limited to purely resistive boundaries, 
but Eqs. (16) and (18) can be generalized to various 
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boundary conditions defined by a reflection coefficient, 
R(uo). This is done by replacing the products like rH{t) 
by the convolution product (r * H){t), where r(t) is the 
inverse FT of R(u>). 

- for the case under study we notice that the convolu- 
tion product of n times function f(t) is (rri) n S(t—2n£/c). 

- if C (respectively (g) is unity, the reflection coefficient 
r (respectively re) vanishes, as well as fit): the first term 
of the Green function is the Green function of an infinite 
medium, the first two terms correspond to a semi-infinite 
medium, etc... As it will be seen in the next section, no 
modes can be found for these cases, because no reflections 
exist, either rj or rji tending to infinity. 

- finally, multiplying both members of Eq. (15) by the 
factor [1 — F(u>)], and taking the inverse FT, it appears 
that a closed-form exists in the time domain, which is 
the basis for the study of the Hclmholtz motion of bowed 
string instruments (see e.g. Woodhouse 18 ). It is a re- 
currence relationship: 

d t g(x,t) ~rred t g(x,t-2£/c) = g p (x,t) c/2. 



ii) if | CI > 1 and \(e\ < 1 (mixed case with either 
complex 77 or complex rji : either /1 or \m is unity) : the real 
part corresponds to a problem with different conditions 
(Neumann and Dirichlet) at x = at x = I. The real 
part of eigenfrequencies is an odd harmonic of c/At. 

iii) if |e| < 1 and \Q\ < 1 (complex r\ and r\t- M = 
jit = 1): the real part corresponds to the values for pure 
Dirichlet conditions (zero £ and Q). 

Except for case ii), a purely imaginary eigenfrequency 
exists for n = {fi + ^e)/2. 

B. Calculation of the residues 

In all cases, the Taylor expansion of the function sinh 
in Eq. (12) at the first order of the quantity (w-w„) can 
be determined. The result is: 

sinh [iuol/c + 17 + ~ i(-l) n (w - u n )l/c. (22) 
We get for lo close to the pole co n : 



V. EXPANSION IN RESONANCE MODES USING 
THE INVERSE FT 

Putting expression (12) of the frequency domain in 
Eq. (5) leads to the modal expansion of the time domain 
expression. The tool is the residue calculus. If all poles 
of expression (12) are simple and located on or above the 
real axis, the following equation can be used: 



g(x,t)=iY, if t > and if t < 0, 



(19) 



where S is the sum of the residues of G(x, lo) exp(iuit) 
(see e.g. Morse and Ingard 8 p 17, changing i to —i). 



G(x, Lo) 



C 2 fn(x)fn(x ) 



LO n £ (lo - LO n ) 

f n (x) = cosh(?7 + iuj n x/c) 



(23) 
(24) 



or f n (x) = (— 1)" cosh(iui n (£ — x)/c + rji). The residue 
corresponding to the pole lo — 0, remains to be calcu- 
lated. For small lo, 

„, . c cosh n cosh % c 1 , . 

G(x,lo) = —i £ = r . 25 

ilo smh(r) + r)i) ilo C,- 1 + Q 

Using Eq. (19), the inverse FT of G(x,oo) is obtained: 



A. Calculation of the poles 

Zeros of function sinh satisfy: 

ito n = [—rj — rje + iti-k] c/£, (20) 

where n is an integer. In order for the poles to be above 
the real axis, the condition is rj r + rie r > 0. It is equivalent 
to the condition previously obtained for the successive 
reflections expansion: \rre\ < 1. Using definition (10), 
Eq. (20) is rewritten as: 

w n = [n - (ji. + m)/2] irc/£ + i(rj r + rj lr )c/£. (21) 

As already remarked by several authors, the imaginary 
part of the complex frequency is independent of n, and 
the real part is independent of the dissipation. Depend- 
ing on the values of Q and Q, different cases must be 
distinguished: 

i) if |C| > 1 and |^| > 1 (real i] and rji): the real 
part of the frequency corresponds to the values for pure 
Neumann conditions (infinite Q and Q). 



g(x,t) = Hitfl £ mi^l e ^ + ^) . 

(26) 

Some comments can be made: 

- the formula is valid for all aforementioned cases; 

- the mode shapes f n (x) are complex- valued functions 
of the space variable, meaning that the shape is vary- 
ing with time (for a discussion on complex modes, see 
e.g.Ref.19). The question of their orthogonality will be 
discussed in section VI. Notice that functions f n (x) do 
not fulfill the same boundary conditions than G(x,lo): 
the boundary conditions are (7) and (8), but where lo is 
replaced by Lo n ; 

- the last term in Eq. (26) is a constant mode. If one 
of the impedances £ or Q is zero, it disappears, as it is 
intuitive, in order to satisfy a Dirichlet condition. It is 
a trivial solution of the wave equation and the boundary 
conditions, and can be compared to the DC component 
of a periodic signal. When both ( and Q tend to infinity, 
the boundaries tend to Neumann boundaries, and the 
combination of the non oscillatory mode of frequency u)q 
and the constant mode results in a uniform (i.e. constant 



4 J. Acoust. Soc. Am. 



Kergomard: Resonance modes in a 1 D medium 



in space) mode which increases linearly with time. The 
calculation is done as follows: if r\ and rj£ tend to zero, ujq 
tends to zero, and the factor exp(iuJot) can be written as: 
1 + iu>ot. The zeroth order term is equal to the opposite 
of the constant mode, and only the linear term remains. 
The result is H(t)tc 2 /£, and its FT is -c 2 /£oj 2 . This 
mode is the classical uniform mode existing for instance 
in 3D cavities with rigid walls: curiously it exists in the 
standard textbooks (see e.g. Ref.8, page 571), but the 
time domain expression is not given. This mode is similar 
to the well known planar guided mode, existing in ducts 
with rigid walls, whatever the geometrical shape. 

- the imaginary part of the complex frequencies be- 
ing independent of n, the decay is identical for all non 
constant modes; 

- for the above-considered case i), we notice that w_„ = 
— oj* and /_„(x) = f„(x), and, more generally: 

iw„ = (iu n )* ; U(x) = (-irf*(x) , (27) 
where v = —n + fi + fig is an integer. (28) 

As a consequence, the solution g(x,t) is real. It could be 
possible to transform the sum by adding the two oscillat- 
ing terms corresponding to n and v, when n ^ v, as it is 
usually done for non dissipative boundaries. Nevertheless 
it appears that the formulas become intricate. 

- Eqs. (23) and (25) lead directly to another form of 
the FT of result (26), written as a series: 



G(x,oj) 



<?_ fn(x)fn(x ) 
£ ^ UJ n {0J - L0 n ) 



(29) 



An example of comparison of the successive reflections 
method and modal expansion is shown in figure 2(a), for 
conditions close to Neumann and Dirichlet. We notice 
that it is satisfactory. The Gibbs phenomenon appears, 
because of the truncated series of modes, ensuring the 
correct accordance between the two methods. Moreover 
this accordance confirms the existence of the constant 
mode. Decreasing of the maxima is due to the dissipa- 
tion at the boundaries: for pure Neumann and Dirichlet 
conditions, the shape would be similar, but perfectly pe- 
riodical. 
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FIG. 2. Normalized Green function as a function of time: 
comparison between the successive reflections method (dot- 
ted line) and modal expansion (solid line, 102 modes, i.e. 
maximum n = 50). Locations of the source and receiver are 
xo/l — 0.15 and xjl = 0.76, respectively, 
(a) passive boundaries f = 4. ; Q = 0.1 (constant mode= 
0.097). (b) active boundary ( = C,i = —4. (constant mode= 
-2). 



and similarly for Eqs. (2) and (3). It is equivalent to 
use an appropriate Laplace Transform. Going in the fre- 
quency domain leads to Eqs. (6 to 8), where G(x,u>) 
is replaced by G(x, to) and iu> by (iu> + rj), and a sim- 
ilar result for Eq. (12). The analysis of both succes- 
sive reflections and poles and residues leads to the re- 
sult g(x,t) — g(x, t) exp(— rjl), where g(x,t) is given 
by Eqs. (17) and (26), respectively, and the proof is 
achieved. We do not repeat here the complete procedure. 
We notice that for the case r\ r + r/£ r = 0, one boundary is 
active and the other one is passive: eigenfrequencies uj„ 
are real while modes are complex. Figure 2(b) shows an 
example of result. 



C. The case of active boundaries 

What happens when the combination of boundaries is 
active, i.e. when \rri\ > 1, or r\ r + r\i r < (at least 
one of the impedances £ or Q is negative)? It is possible 
to prove that Eqs. (17) and (26) remain valid for active 
boundary conditions, as explained hereafter. The real 
part of iuj n being independent of n, this suggests to use a 
new function ~g(x, t) = g(x, t) exp(-rjt), where rj > — r/ r — 
r]i r > 0, which can be substituted in the initial problem 
in order for the poles to be located again on or above the 
real axis. Eq. (1) becomes 

d 2 xx g(x, t) - c- 2 [d t + g(x, t) = -8(x - x )e-^6(t) 



VI. METHOD OF SEPARATION OF VARIABLES 

A. Second order homogeneous equation with 
initial conditions 

1. Derivation of the modes 

Oliveto and Santini, 2 and Guyader 3 have treated a 
particular case of the problem (zero ( , large Q) us- 
ing the method of separation of variables. He gets 
non orthogonal modes for the common scalar product 
Jq fn{x)f m (x)dx. We will see that the method is valid 
whatever the values of the two boundary conditions, and 
how the derivation can be simplified. 
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We are searching for solutions p(x, t) of homogeneous 
equation (1) (without second member), with boundary 
conditions (2) and (3), and with given initial conditions. 
Assuming that the general solution is a superposition of 
solutions with separate variables, the solutions with sep- 
arate variables are written in the following form: 



p(x,t) = f(x)h(t); 
h(t) = B+e iut + B-e~ ilvt ; 
f(x) — cosh(iujx/c + f). 



(30) 
(31) 
(32) 



Decomposition (30) differs from the ordinary FT, be- 
cause a priori lo is a complex quantity, depending on 
the boundary conditions. Considering first the solution 



B+ 



this leads to : 



(lo sinh ip = lo cosh <p; (33) 
u)Q smh(ioj£/c + lo) = — lo cosh(iui£/c + p). (34) 

lo = is a solution, corresponding to the constant 
mode. The other modes are given by Eq. (33): smh(p> — 
i]) = 0, thus: 



f(x) = cosh(iLox/c + rj). 



(35) 



Actually there is a sign ± in the right-hand side mem- 
ber of Eq. ( 35), but it is without importance, because 
it can be included in the coefficient B + of the solution. 
The eigenvalues equation is deduced from Eqs. (33) and 
(34), as follows: 



sinh(ikj^/c + rj + 7]e) = 0, 



(36) 



the solutions being given by (20). The solution in time 
g- e -iu>t ^ogg not i eac j £ Q new som ti ons for f(x), there- 
fore, assuming the solutions form a basis of solutions (this 
is discussed in section VI B) , the general solution of a 
problem with initial conditions can be written as: 



p( X ,t) =J2 An/nix^* + A, 



(37) 



where ui n and f n {%) are given by Eqs. (20) and (24), 
respectively, and the coefficients A n and A depend on 
the initial conditions, and can be determined using the 
orthogonality relation of the modes. A is the coefficient 
of the constant mode. 

2. Orthogonality relationship between the modes: 
first approach 

In order to derive an orthogonality relationship be- 
tween the modes the common product is first calculated: 



A n m / fn (*^)/m {x)dx. 




(38) 



Because f n (x) = (—l)^f*(x), the calculation of the 
quantities defined in Eq. (38) for all values of the index n 



is equivalent to the calculation of the quantities defined 
when replacing f m (x) by its conjugate. Writing 



dx 2 



dx 



dx 2 



and using Eq. (24), the following result is obtained: 



- ul) 



A, 



l(LJ m -LO n ) 



/n(0)/ m (0) , f n (£)f m (t) 



C 



For Lo m ^ Lo n , because Lo m + Lo n ^ 0, the expression 
of A nm is deduced. For Lo m = u> n , the calculation is 
straightforward. The general formula is found to be: 



CI 



LO m + LO r , 



/n(0)/m(0) fn(l)fm{£) 

C St 



1 



where S nm is the Kronecker symbol, or: 



(39) 



c smh2r,+ (-1)"+- sinh 2 m 1 

2 l (LO m +LO n ) + 2 ^nrn- (4U) 

Modes are found to be non orthogonal for the product 
defined by (38), but, as shown by Guyader, 3 it is pos- 
sible to solve the problem from the knowledge of initial 
conditions. When dissipation tends to zero (r] r and r\i T 
tend to zero), the first term does not vanish, tending 
to {—lY\£8 n ^. This is due to the choice of considering 
separately the modes u> n and lo v . 

Otherwise formula (39) remains valid when one of the 
modes is the constant mode f(x) = 1, and the other one 
a non constant mode: 



/n(0) , fn(£) 



A n = / f n (x)dx = 

Jo l ^n 

= —(sinh 7?+ (-l) n sinh^). (41) 

lLO n 

Finally the product of the constant mode by itself is £ . 



3. Solution with respect to initial conditions 

According to Eq. (37), the initial conditions are: 

p(x, 0) = A n cosh(iLo n x/c + 77) + A; (42) 

n 

d t p(x,Q) A n iw n cosh(iw n x/c + rj). (43) 

n 

Using Eq. (40) for a non constant mode m, the fol- 
lowing results are obtained: 

rt 

/ p(x,0)f m (x)dx =^TAnAnm + AA m ; (44) 

rt 

I d t p(x,0)f m (x)dx =\ (45) 
Jo „ 
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Multiplying Eq. (44) by iui m , then adding Eq. (45), 
leads to: 



iP(x, 0) + d t p(x, 0)] f m (x)dx 



uj n )A, 



/ m (0)p(0,0) + / m (*)p(*,0) 



C 



6 



+ iAu m A m 
+ iA m £uj n (A6) 
+ iA m £u m (A7) 



Notation f° r tnc scr i cs in Eq. (46) indicates that 
it involves the constant mode. As noticed by Guyader, 3 
this series is related to the initial conditions at the two 
ends x = and x = £. Thus for a non constant mode: 

A n £iuj n = / [iu n p(x, 0) + d t p(x, 0)] f n (x)dx 
Jo 

+cp(0, 0) sinhr/ + cp(£, 0)e„ sinh%. (48) 

The following property is deduced from Eq. (27): 
A v f v {x) = A* l f*(x), thus p(x,t) is real. Calculating 
Jo dtp(x, 0)dx, we similarly get coefficient A: 



A = 



c 1 J o d t p(x,0)dx + p(0,0) tanhrj + p(£,0) tanhrje 

tanh r\ + tanh rje 

(49) 

What is the condition for which this coefficient van- 
ishes? If for instance at x = 0, (is zero, rj is infinite, 
and, according to the boundary condition, p(0, 0) van- 
ishes, thus A vanishes too. This confirms the remark 
concerning result (26). 

Using the initial conditions for the Green function found 
in (17), it is possible to check result (26), but this will be 
done hereafter using the equation with source. 

B. First order system of equations, orthogonality 
and completeness of the modes 

1. Introduction 

In this section we will prove that the modes form a 
Riesz basis in the space of solutions of a closely related 
problem, and give the expression of a scalar product mak- 
ing the modes orthogonal. As an introduction we show 
that a modified scalar product leads to the orthogonality 
of modes, except the constant one. For vibrating systems, 
the product defined by (38) corresponds to the product 
with respect to the mass, a complement being the cal- 
culation of the product related to the stiffness (see e.g. 
Meirovitch 20 ): 

A ' nm = J i fn ^iL fm( - x ' ,dx ' 



By integrating by parts, and using Eq. (39), this prod- 
uct, for n 7^ m, is found to be equal to: 

2 

A' =\f (x)d f (x)f + ^-A - UJnUJm A 

Therefore the modes become orthogonal if we define a 
new product, as follows: 



o 



$xPn &xPrt 



1 



dtPndtPn 



dx = S nm £ul/c 2 



t=0 



(50) 

where p n — p n (x,t) — f n {x) exp(iw n t) and similarly 
for index m. We remark that the modes p n and p v = 
(— l)^Pn are orthogonal for this product. For the calcula- 
tion of the solution from initial conditions, using Eq. (37) 
at t = 0, the following result is obtained: 



-^fn(x)d x p(x,0) - ^-f n (x)dtp(x,0) 



dx — A n <. 



M. 
c z 

(51) 

As a consequence, the initial conditions need to be 
written by using the derivatives of the function p{x, t) 
with respect to abscissa and time, respectively. Result 
(48) can be checked by integrating by parts the first term 
of the integral. Nevertheless, the product (50) is not use- 
ful for the constant mode, and the first method needs to 
be used (see subsection VI A3). Moreover this derivation 
does not prove that the product is a scalar product, and 
that the modes form a basis of the space of solutions of 
the problem. This will be done hereafter. 

2. Riesz basis of the modes 

Several works have been done by mathematicians 
concerning spectral operators when boundary condi- 
tions are not simple conditions like Neumann or Dirich- 
let conditions. We quote the work by Russell, 21 
Majda, 22 Lagnese, 23 Banks et al, 24 Darmawijoyo and Van 
Horssen, 7 Cox and Zuazua. 6 Rideau 4 has treated the ID 
case with a (unique) resistive termination, giving explic- 
itly a scalar product (see also Ref.5). We generalize his 
calculation using a similar method, by considering the 
wave equation with source in the following form: 



d t tp(x, t) = Atp(x, t) + <j> s (x, t), 



(52) 



where tp(x,t) — (p,v) T , p and v/(pc) being the acoustic 
pressure and velocity, respectively. Operator A is: 



A = 



-cd x 
-cd x 



(53) 



and boundary conditions are written as: 

p(0,t) = -£v(0,t) and p(£,t) = Qev{£,t) Vi. (54) 
The family of eigenelements of A are found to satisfy: 
KPn(x) = -cd x v n (x) ; X n v n (x) = -cd x p n (x), (55) 



J. Acoust. Soc. Am. 



Kergomard: Resonance modes in a 1 D medium 7 



thus 



Pn(x) 

v n (x) 



cosh(A„x/c + rf) 
— sinh(A„x/c + rf) 

r]t + inn) c/£ - 



-n 



(56) 
(57) 



(see Eq. (20)). p n {x) = f n (x) and A„ are identical to the 
cigcnfunctions and eigenvalues found before. Neverthe- 
less the constant mode is eliminated (except for the very 
particular case r\ = —Tfe), because the boundary condi- 
tions are slightly different: Eqs. (2) and (3) are obtained 
by deriving Eqs. (54) with respect to t. In Eq. (56) the 
argument of the hyperbolic functions can be written as: 

X n x/c + r] = a(x) + i/3 n (x); (58) 
a(x) = — rj£ r x/£ + i] r (1 — x/£) ; (59) 
(3 n {x) = 7T [-fj, e x/£ + fi (1 - x/l)] /2 + mra/<60) 

(see definitions (10) and (11)). Denoting ip"(x) = 
(p n {x), v n (x)) T , we show in Appendix A that the family 
of elements ip%(x) is a Riesz basis, i.e. a complete basis 
of elements, which become orthogonal for the following 
scalar product: 



<xjj lV > a H = / <p T *G 2a (x)ipdx., 
Jo 



where G a (x) 



cosha(x) sinha(x) 
sinha(x) cosha(x) 



(61) 



(62) 



For a given vector tp = (p 7 v) T , the value of the scalar 
product with eigenvector is found to be: 



< 



fp,^n >h= / \p(x,t)p n (x)-v(x,t)v n (x)]dx. (63) 
Jo 



This is in accordance with the product (50). A di- 
rect application of this result is the solution of Eq. (52) 
with initial conditions ip(x,0) — (p(x, 0), v(x, 0)) T ■ 
The modal decomposition is uniquely determined as 
ip(x, t) — J2 h n (t)tp"(x) in the energy space H, and leads 

n 

to the following family of decoupled ordinary differential 
equations: 

£ [dth n - \ n h n ] = < 4> a (x, t), c >h; (64) 

MO) - <^(x,0),C>h- (65) 

This result can be first applied to the calculation done 
in section VIA. In order to find a solution x(x,i) of 
the second order equation without source, we denote 
ip(x,t) = (dtx(x, t), — cd x x{x, t)) T , and obtain by inte- 
grating ip(x,t) with respect to t : 



X (x,0) =J2 K 1 hn(0)Pn(x)+A. 



(66) 



Using Eqs. (65), (63), and replacing x( x ,t) by p(x,t), 
Pn{x) by f n (x), \ n by iw n , and A" 1 /i„(0) by A n , for- 
mula (51) is checked. Notice that coefficient A cannot be 
directly determined with this method. 



3. Example of the Green function 

Similarly, the Green function can be calculated by us- 
ing the previous result. In order for the unknown function 
to satisfy the boundary conditions (54), or (2) and (3), 
it is convenient to define the following vectors: 

«*■«-( -28$) 

(67) 

The first row of Eq. (52) is Eq. (1), while the second one 
comes from the definition of vector xjj. Using Eq. (64), the 
solution is found to be: ip(x,t) V'n ( x )h n {t) 7 where 

n 

d t h n - X n h n = c 2 £~ 1 p n (x Q )S(t). (68) 

The initial conditions for the Green function imply 
tp(x,t) = for t < 0, therefore h n (t) = for t < . 
Thus the solution of Eq. (68) is: 

K{t) = A n e Xnt H(t); A n = J^p^). 

As a consequence, 



d t g(x, t) 
-cd x g(x,t) 



Pn(x) 

v n (x) 



e A "*iJ(t). (69) 



Integrating the first row with respect to time leads to: 
g(x, t) - H(t) [A nPn (x)e x - t + A{x)\ . (70) 

n 

Dcrivating this expression with respect to x and using 
the second row of (69) leads to d x A(x) = , thus A is 
a constant, as expected. In order to deduce the value of 
this constant, we need the following result: 

d t g(x,0) =p(x,0) =c 2 5(x-x Q ). (71) 

It is obtained by derivating the first row of (69) with re- 
spect to time, and the second row of (69) with respect to 
abscissa, leading to p(x, 0)6 (t) = c 2 6(x — x )6(t) (remind 
that d t [F{t)H{t)\ = H(t)d t F(t) + F(0)6(t) ). The end of 
the calculation is done in section VI A3, giving Eq. (49), 
by replacing p(x,0) by g{x,0) and taking into account 
that g(0, 0) = g(£, 0) = 0. We notice that the calculation 
is valid for both passive and active boundaries. 

VII. EIGENMODES EXPANSION IN THE 
FREQUENCY DOMAIN: BIORTHOGONALITY 

Frequency domain approach is very popular in acous- 
tics (see e.g. Ref.8 ), and leads to the use of biorthog- 
onality (see e.g. Ref.17, p. 884) of modes, except when 
the boundary impedances are imaginary, corresponding 
to non dissipative boundaries: for that case, modes are 
orthogonal, and the laplacian operator is self-adjoint. In 
this section we limit the discussion to the Green func- 
tion calculation, and use successively the two above-used 
approaches : the second order equation, then the sys- 
tem of two first order equations, ignoring the constant 
mode. Because we are now in the Fourier domain, equa- 
tions are ordinary differential equations, biorthogonality 
theory ensuring the completeness of the modes family. 
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A. Solution of the second order equation 

1. Modal expansion 

In order to calculate the inverse FT of G(x, uj), another 
solution is possible: the expansion of G{x,uj) in eigen- 
modes. This is done for a particular case by Filippi 1 (p. 
58 : this author considers another type of excitation in- 
stead of the Dirac function, thus uses the Laplace Trans- 
form instead of the FT; notice that the constant mode is 
missing in this work) . We will see how this method leads 
to the same poles and residues that the direct method 
using the closed-form expression (12). We are searching 
for the following expansion : 



G(x,u) G n (x,uj), 



(72) 



where the eigenmodcs G n (x,uj) are solutions of: 

[5L + ^H/c 2 ]G„(x,c)=0, (73) 

and satisfy the boundary conditions (7) and (8). The 
key point is that eigenmodes G n (x,oj) and eigenfrequen- 
cies 6 n (w) depend on frequency w : this is due to the 
boundary conditions, which are of Robin type. Solutions 
of Eqs. (73) can be written as follows: 

G n (x,w) — cosh(id n (uj)x/c + <p n (u)) (74) 

where 9 n (u) and <p n (w) are given by the boundary con- 
ditions. Thus they satisfy: 

6 n (uj) tanh<£„(u;) — u/(; (75) 
O n {uj) tanh(i6 n {uj)e/c + tp„{oj)) = -oj/Ce- (76) 

Eliminating quantity ip n {oj), the eigenvalues are found 
to satisfy the following equation: 



tanh(i#„(w)£/c) 



9 n {uj) + 





'1 1 " 


= —uj 





(77) 

When 9 n (uj) and uj are not simultaneously zero, this 
equation can be rewritten as: 



a 2i8 n e/c _ 



9 n ( + uj 



9n(e - u 



(78) 



Calculation of all solutions of this equation is not nec- 
essary, only one of them being useful in the following. 
Operator D = d 2 x is formally equal to its adjoint D, but 
the boundary conditions are different (conditions for D 
are complex conjugate of conditions for D). Modes of 
D are the complex conjugate of modes G n (x,u>) (they 
are equal to modes G n (x,uj) only if ( and Q are imagi- 
nary, because of the factor i in boundary conditions (7) 
and (8)). Thus in general operator D is not self-adjoint, 
and eigenmodes G„(x,uj) and G n (x,u>) — G* n {x,uj) are 



biorthogonal (see Ref. 17 ). The scalar product of modes 
G n (x,uj) with G m (x,uj)is simply given by: 



/' 

Jo 



G m {x,uj)G n {x,uj)dx = T n S„ 



where 



r„ 



l + c 



sinh 2(i9 n £/c+ <f n ) — sinh 2ip n 
2i£9 n 



(79) 



(80) 



Therefore modes G n (x,w) are orthogonal (for the 
product (79)) and fulfill the same boundary conditions as 
G(x,uj), contrary to resonance modes f n (x) in Eq. (29). 
Finally the solution of Eq. (6) can be written as follows: 



G(x 1 lu)=c 2 J2 



G n (x,uj)G n (xo,uj) 



(81) 



2. Calculation of poles and residues 

In order to calculate the inverse FT, the residue cal- 
culus will be used again. The only terms of the series 
contributing to poles verify: 



8 n (uj) = ±uj. 



(82) 



Looking at Eq. (78), it can be seen that these two solu- 
tions lead to the same equation for u. Rewriting Eq. (78) 
by using Eqs. (75) and (76), the resonance modes fre- 
quencies are found to be solutions of Eq. (20). Solutions 
uj p of this equation are the non zero poles of the integral in 
the inverse FT. Nevertheless the pole uj = exists again, 
because the zero value satisfies Eq. (82), the eigenvalue 
8 n (uj) = satisfying Eq. (77). 

It remains to calculate the residues. Starting with the 
poles uj p ^ 0, we need to select in the series (81) the terms 
involving poles. For a given u> p , there are two terms. 
However it appears that modes corresponding to 9 n and 
— 9 n are identical. As a consequence, only one term of 
the series contributes to the inverse FT: it will be denoted 
9 p (u). The corresponding residue is found by expanding 
Eq. (81) for uj close to ui p , as follows: 



G(x, uj) 



G p (x,ujp)G p (x ,uj p ) 



r p (2u p ){u, p -w){i-[£e p {w)] u=up y 



A similar expression can be found in Filippi, 1 which 
points out that Morse and Ingard 8 (p. 559) forgot the 
derivative. The same error is found in Morse and 
Feshbach 17 (p. 1347), with another error in the deriva- 
tion of Eq. (77): these authors treated the problem of a 
string with one non-rigid (and resistive) support. 

Actually the derivative of 9 p (uj), denoted 9' p (uj) can be 
calculated analytically, as follows. Taking the derivative 
of Eq. (78) with respect to w, or more conveniently, tak- 
ing the logarithmic derivative of Eqs. (75) and (76), the 
following results are obtained: 
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where 



1 0' P _ 2<p' P _ 2(i9' p £/c+cp p ) 
Lu p Op sinh2<p p sinh 2(iQ p £/c + ip p ) 

Thus, eliminating the derivative ip' p , writing O p = lu p 
and using Eq. (80), it is found after some algebra: 

1 - e' p = £/2T p . (83) 
Finally, for to close to co p ,: 

G{x,u) = _* Gp(*,»p)G p (xo,» p ) (84) 

I LUp(UJ - LUpj 

which is in accordance with Eq. (23). 

Otherwise, for u> close to 0, the solution 9(io) which 
is close to 0, solution of Eq. (82), satisfies the following 
equation, deduced from (77) : 

[i - e 2 /3 + o(e 4 )] [e 2 + fi 2 /(CO)] = in [c 1 + Q 1 } 

where = 81 /c and £1 = ujI/c. Therefore 6 2 is of 
order 0, and 

9 2 = itucr 1 [C 1 + Q 1 ] + 0(lu 2 ). (85) 

Using Eq. (81), the residue for the pole u> = is ob- 
tained, and Eq. (25) is confirmed. We conclude that 
the method of the expansion in orthogonal modes in the 
Fourier domain leads to the same result (Eq. (26)) than 
the "direct" method, but the derivation is more delicate. 

B. System of two first-order equations 

We consider now the FT of (52): 

iuitp(x, u>) = Atp(x, u>) + 4> s {x, u>), (86) 

where &(x, to) = (P{x, lu), V(x, lu)) t . An interest of the 
system is that the boundary conditions are independent 
of frequency: 

P(0, lu) = -(V(0, lu) ; P(t, lu) = QV(£, u). (87) 

Eigenvalues A„ and eigenvectors of operator A arc 
already known (see Eqs. (56) and (57)). The appendix 
shows that the adjoint operator of A, is A = —A, and 
gives the boundary conditions for it. This formulation 
differs slightly from the work, 10 ' 11 considering a different 
operator, but the principle is identical: we notice that 
these authors treat the problem for more general opera- 
tors and boundary conditions. Eigenvalues and eigenvec- 
tors of A are solutions of: 

= (v^ for x = p^l = -Qu^ for x = I. 

Thus, the adjoint eigenvalue problem to be solved is 
the same as the direct one, by replacing c by — c, r\ by -r\ 
and rje by —rig. The eigenelements are thus found to be: 

tut _ ( P^(x) \ _ ( cosh(A m a;/c + r?) \ , . 

^ m { v^(x) ) \ sinh&/c + rj) ) ( ™ } 
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A m = - (r/ + r\i + imn) c/£ = A_ TO (89) 

Comparing with the family tp" (Eq. (56)), there is 
a difference in sign for the second row: we notice that 
Rideau made an error in the biorthogonal family. By 
construction, the biorthogonality relationship is ensured: 

(A„-A^*) <C^>=0. (90) 

Using Eq. (27), we remark that A„ = A m * implies m — 
—v, as defined in Eq. (28). Therefore 

<^M>= f^ T *^dx=(-lTt5 rn ,_ v . (91) 
Jo 

This latter relation enables to perform a modal de- 
composition on the (ip") family: but, contrarily to 
standard cases, the n-th coefficient is not given by the 
scalar product with but by the scalar product with 
i^-i/ ! ( U P to the normalization coefficient (—1)^1). 
Notice that if rj and r/g are both real, A m = A*„, and 
(91) is obvious from the expressions of eigenvalues 
and eigenvectors (for this case, A„ = A m implies 
n = m). For the general case, the scalar product can 
be written: < ip^ >= J Q cos [(3 n (x) + (3- m (x)] dx = 

(—l) fJ, J cos[(i' + m)TTx/£]dx. Comparison with 

Eq. (A4) exhibits the difference between the two meth- 
ods. 

It remains to apply orthogonality to Eq. (86). We choose 
the case of the Green function (Eqs. (67)), with the 
following result: 

GM = - c l y fjMIM. (92) 

The calculation is easy, because u> n docs not depend on 
frequency. Comparison with Eq. (29) exhibits a differ- 
ence in the denominator, i.e. a factor lu instead of aj n , 
and, of course, the absence of constant mode. When re- 
turning to the time domain, all the terms corresponding 
to to = to n are identical, and a constant mode is found 
for the pole u> = 0, but again it is not possible to de- 
duce it from orthogonality relations, as in section VI B. 
Nevertheless, because of the independence of the bound- 
ary conditions with respect to frequency, the calculation 
of the residues is much easier than for the second-order 
equation. For the same reason, the calculation in the time 
domain would be possible with the same modal decom- 
position, and this is a major difference with the methods 
based upon the second-order equation. 

VIII. CONCLUSION 

The simple problem we have studied, which can be 
regarded in particular as a radiation problem, exhibits 
interesting properties for the resonance modes: they are 
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complex- valued, and non orthogonal for the simple prod- 
uct (38) because of the bounded character of the con- 
sidered medium, but except the constant mode, they are 
orthogonal for a product modified in a proper way, and 
are a basis for the space of solutions. Second order equa- 
tions allow to find the constant mode, while first order 
systems of equations allow a more direct formulation of 
boundary impedances. 

Thanks to the simplicity of the problem, the analytical 
treatment is possible with several methods, elucidating 
the relationship between them, which can be useful for 
more intricate problems (e.g. when damping is added 
to propagation, or when boundary impedances involve a 
mass or a spring). No approximations arc needed, the 
results are valid whatever the value of the terminal resis- 
tances. Active boundaries can also be considered, thanks 
to a change in functions. We notice that an advantage 
of the frequency domain calculations is the possibility of 
the treatment of an arbitrary dependence of the bound- 
ary conditions. For a dependence rj(u)) and %(w), Eq. 

(29) remains valid by replacing £ by I — ic(r/' p + r)' lp ) , 

where rj v = (dr//duj) u)=LUp , and similarly for r/g. This can 
be shown by generalyzing Eq. (22), or, with some alge- 
bra, using the modal expansion. 

Finally, considering the problem of a stratified medium 
(see section II), it could be deduced in the field outside 
of the interval [0, £] . When terminations are passive, a 
result is that modes tend to infinity when x tends to 
±oo. An interesting study has been done in Ref.14, using 
biorthogonality and explaining the relation between the 
energy outside the interval and the terms responsible of 
non orthogonality in equation (39). 
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APPENDIX A: PROOF OF THE COMPLETENESS OF 
THE EIGENELEMENTS OF OPERATOR A 

The operator A, defined by Eq.(53) is a differential op- 
erator defined on the energy space H = L 2 (0,£)xL 2 (0,£), 
it has a compact resolvent (cf Ref.5, p. 191). Using the or- 
dinary scalar product < xp, if >= J Q [pq* + vw*] dx, be- 
tween ip(x, t) = (p, v) T and <f(x, t) — (q, w) T , the follow- 
ing result is obtained: 

< Atp, tp > + < ip, Aip > = 

-c [v(q* + <Uw*)] x=t + c [v(q* - (w*)] x= 



o ■ 



It is deduced that the adjoint operator of A is A = A 
(we denote all quantities related to the adjoint problem 
with an overline), and on its domain, the following ad- 



joint boundary conditions must be fulfilled: 

q(0,t) = (*w(0,t) and q(£,t) = -Qw{t,t) Vt 

(here ( — (* , and Q = Q; if £ is infinite, the boundary 
conditions are v(0,t) — 0, and w(0,t) = 0, and similarly 
for boundary x — £). Therefore A is skew- symmetric, 
but not skew- adjoint, because the domains of A and A 
are different, except if both £ or Q are either zero or 
infinite (Dirichlet or Neumann conditions); notice that 
for a skew-adjoint operator, the eigenvalues are purely 
imaginary. In order to find a new scalar product, we 
denote, from Eqs. (58) to (60): 



e a(x) e -a(x) 



i(3„(x) 
-i/3„(x) 



In H the standard scalar product < ip",ipp >h= 
J'q [PnPp + ^n^p] dx does not vanish for n ^ p, except if 
a(x) — 0. If we denote ipn(x) the functions corresponding 
to the latter case, it is possible to construct a new scalar 
product ensuring orthogonality, in a similar way Rideau 4 
did. From Eq. (62), the following hyperbolic rotation is 
obtained: 



^W = G„(x)C(4 

We will now prove that the new product 



<V >H= / V 
JO 



(Al) 



M a (x) xp dx 



<ip,(p>%=<G a ip,G, 

(A2) 

where M a (x) = G^G a , leads to the orthogonality of 
the modes. M a (x) is found to be equal to G^a- It is 
symmetrical and positive definite because 



(NCr = ( (p,« 



HQ 2 



cosh [2a(x) 



can be rewritten as: 



H 2 ) + 2sinh [2a(x)} $e(pv*) 



dx 



(II 



\ a H f 



1 

2./,, 



\p + v\ 



e -2a(x) \ p _ v f 



dx. 



Moreover a(x) is a function varying monotonously 
from r\ r to — r\i r when x increases from to £, and the 
following bounds can be found for HVHIh 1 



C a \\M\h < UWl < C a \\M\h 



(A3) 



where c a — e~^ and C a = e v , with rj — sup [\i] r \ , \r/e r \}. 

Therefore the modes ip„ are orthogonal for the new 
scalar product < xp,tp > c f { . First recall that (ipn) n i s 
the family of eigenvectors of a classically skew-adjoint 
operator with compact resolvent, it is thus complete in 
H. Now, thanks to (A3), the two norms are equivalent 
on H, and the hyperbolic rotation shows that (ipn) n anc ^ 
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(ipn)n s P an the same subspace, namely the whole of H . 
This proves the completeness of (■*/>")„ in H. 
The calculation leads to the simple result: 

< ^n^m>H = <Vn^m>H= / COS [/?„ (x) - (3 m (x)] dx 

Jo 



i 

Jo 



cos [{n — m)TTx/£] dx = £5 nm . 



(A4) 
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